Exciton effects in a scaling theory of intermediate valence and Kondo systems 
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An interplay of the Kondo scattering and exciton effects {d — f Coulomb interaction) in the 
intermediate valence systems and Kondo lattices is demonstrated to lead to an essential change of 
the scaling behavior in comparison with the standard Anderson model. In particular, a marginal 
regime can occur where characteristic fluctuation rate is proportional to flow cutoff parameter. In 
this regime the "Kondo temperature" itself is strongly temperature dependent which may give a 
key to the interpretation of controversial experimental data for heavy fermion and related systems. 
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There is an interesting class of rare-earth compounds 
such as Ce4Bi3Pt3, SmBg, SmS under pressure (the 
"golden" phase), TmSe, YbBi2 which were called ear- 
lier intermediate valence (IV) compounds and now are 
treated as "heavy-fermion (HF) semiconductors", or 
"Kondo insulators" (for a review, see RefsJ., ^). Vari- 
ous names emphasize different peculiarities of these com- 
pounds. As for electron energy spectrum, most of them 
are narrow gap semiconductors with an anomalously 
small energy scale (the gap width) , of the order of tens 
or hundreds of Kelvins (see RefQ and a review of ear- 
lier experiments in Ref. 13). At the same time, they do 
demonstrate intermediate valence of the rare earth ions 
(usually between 2+ and 3-I-) in a number of properties, 
e.g., lattice constants (which are intermediate between 
those for isostructural compounds with di- and trivalent 
ions) , core level spectra (which are a mixture of the spec- 
tra of di- and trivalent ions with comparable weights), 
and many others 

As well as for the HF metals, the origin of this small 
energy scale is a key point to understand anomalous prop- 
erties of the IV compounds. For the HF metals it is com- 
monly accepted now that they are the Kondo lattices, 
which means that this energy scale (the Kondo temper- 
ature Tk) is a width of the Kondo resonance owing to 
spin-dependent scattering of conduction electrons by /- 
electron centers @. As a result of an interplay of the 
Kondo effect and interatomic magnetic interactions, the 
Tk value for a lattice can strongly differ from that for an 
isolated impurity 0, spin fluctuations being of crucial 
importance for the HF behavior. It is very natural to ex- 
pect that these effects are important also for the IV com- 
pounds. At the same time, valence or charge fluctuations 
should be also considered. They are determined in part 
by the Coulomb ( "Falicov-Kimball" ) interaction between 
conduction and localized electrons Q. Taking into ac- 
count these interactions together with the hybridization 
processes it is possible to describe the IV state as a kind 
of exciton condensation [^|^. Note that in the IV regime 
the spinless Falicov-Kimball model with hybridization is 
formally equivalent to the anisotropic Kondo problem, 
different valent states playing the role of pseudospin "up" 
and "down" states [ifl . It is the degeneracy of quantum 



states for a scattering center which is important for the 
formation of the Kondo resonance j^. In the IV case 
the divalent and trivalent states are degenerate by defi- 
nition, so that this analogy is not surprising. Therefore 
it is natural to consider the formation of the Kondo res- 
onance for the IV compounds taking into account both 
spin and charge fluctuations, or, equvalently, both the 
"Kondo" and exciton ("Falicov-Kimball" ) effects. This 
is the aim of the present work. Since there is no clear 
demarcation between the IV and Kondo systems, it will 
be shown that the excitonic effects may be relevant also 
for the latter case. 

To investigate effects of interaction of current carriers 
with local moments we use the Hamiltonian of the asym- 
metric infinite-?/ SU (N) Anderson model with inclusion 
of the Falicov-Kimball interaction (on-site d — f Coulomb 
repulsion G), 

km 

+Efflmhm\+G J2 fLf^m4^,C^m' (1) 
imm' 

where the on-site f — f Coulomb interaction is put to 
infinity, so that doubly occupied states are suppressed, 
fin ~ \'i'm){iO\ are the Hubbard operator (|zm) and \iO) 
are single-occupied and empty states), we neglect for sim- 
plicity the k-dependence of the hybridization V. Note 
that similar calculations can be performed for realistic 
rare-earth ions, including the case of two magnetic con- 
figurations, see Ref. 1% 

Following to Ref. llCi we treat the coherent and inco- 
herent cases. In the first case dispersion in the spectrum 
of /-electrons occurs. For simplicity this is supposed 
to be proportional to the conduction electron spectrum, 
Efk — Ef + Ask, £k ik, A = —1 (the /-band has a hole 
character). In the incoherent regime A = 0, so that the 
/-electrons remain localized. Note that in the presence of 
the energy gap we always deal with the coherent regime. 

The renormalization of the Coulomb parameter G and 
the hybridization V is obtained, similar to Ref. from 
the two-particle Green's function 
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which determines the vertex. We obtain the singular cor- 
rection with the structure 
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where num = (cj^mCkm) is the Fermi function. In the co- 
herent regime, a similar correction occurs from the dis- 
persion of /-states. The correction (jSJ contains a loga- 
rithmic Kondo-like divergence owing to charge fluctua- 
tions, which is cut at Ef (the latter quantity plays the 
role of the external field in the equivalent anisotropic 
Kondo model). Unlike the renormalization of Ef, the 
renormalizations of V and G do not contain the degen- 
eracy factor of N. 

The renormalization of Ef owing to spin- flip processes 
is obtained in the second order in hybridization (cf. 

Iki8j) 
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We have taken into account in Eq.Q the Hartree renor- 
malization of /-level energy, which occurs in the coherent 
case, Ef ^ Ef + On, n being the concentration of con- 
duction electrons; we put in numerical calculations n = 1, 
which corresponds to the Kondo regime. 

To derive the scaling equations for the effective model 
parameters we use the poor- man scaling approach [T3 |. 
Picking out in the integrals with the Fermi functions lpj)l. 

the contributions from the energy layer C < E < 
C + SC near the Fermi level Ep — and replacing Ef 
Ef{C), V ViC),G G{C) we obtain (cf. Ref. Efl) 
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where p is the bare conduction-electron density of states 
at Ep. Earlier i| we have considered the exciton effects 
with neglecting spin fluctuations. We will see that the 
renormalization results in new essential effects. 
We have from Eqs.® and l(7)) 
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so that G(C) — G{0) in the incoherent regime. We de- 
rive from © in the incoherent and coherent cases, re- 
spectively 

V{G) = ViO)\D/w{C)\P^^"^ (9) 
ViG) = ViO)/[l + G{0)p\n\wiG)/D\] (10) 

where w{C) = G — Ef{G), D is a cutoff parameter of 
the order of bandwidth (we put in numerical calculations 



D = p ^ = 1). Then we have the closed equation for 
w{G). In particular, for the incoherent regime 
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When Ef lies sufhciently below the Fermi level (the 
Kondo regime), the quantity |'u;(C)| can become small 
with decreasing |C|. We can use this condition to define 
the boundary between IV and Kondo cases. Formal def- 
inition of IV systems is the absence of solutions to the 
equation w{G) = which just determines the Kondo res- 
onance (cf. Refs. Ig UT^ . Physically, Kondo lattice has a 
three-peak density of states (two Hubbard bands and the 
Kondo resonance), which is similar to the "doped Mott 
insulator" (note that in the dynamical mean-field theory 
(DMFT) the Hubbard model is reduced to the Ander- 
son impurity model 15]). On the other hand, IV state 
is similar to the phase of strongly correlated metal: the 
Kondo peak as a separate solution is absent. 

For G(0) — 0, the boundary condition for the Kondo 
state is 1^^/(0)1 > T = {N - l)pV^{0). In the opposite 
IV case \w{G)\ remains finite. For G(0) ^ 0, V{G) in- 
creases during the renormalization process, and the ef- 
fective level width r(C) becomes larger, so that the IV 
region becomes more wide. A temperature dependence 
of the energy gap (an increase with decreasing temper- 
ature) in IV compounds was observed experimentally in 
SmBg jlfii, YbBi2 and Cc4BiPt3 2]. According to our 
treatment, the dependence of the effective hybridization 
V{G) is non-monotonous: it passes through a maximum. 

Now we consider in more detail the incoherent case 
which should be realized for diluted systems (the An- 
derson's localization prevents coherence at low tempera- 
tures). To present numerical results (Figs. 1-4) we use the 
variable £^ =\n\D/G\. As follows from ©, in the Kondo 
case the hybridization parameter V{G) decreases practi- 
cally by a jump when we approach the point C = Ef{G), 
see Fig. 2. With further decreasing |C|, a considerable re- 
gion arises where we have to high accuracy G = Ef{G) 
(Fig.l). More exactly, we have dw{G)/dG ~ 0, so that 
we obtain from Hll|) near the maximum of w{G) 
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In this regime V{G) = \G/[{N ~ l)p]\^^^ , and the 
effective s — f exchange parameter is pI{G) = 
pV^{G)/Ef{G) = ~l/iN - 1) = const. 

In a standard consideration, the condition G = Ef{G) 
determines an energy scale for a crossover to the regime 
of a heavy- fermion (Kondo) local Fermi liquid (Ij. The 
"marginal" situation with Ef{G) = C in a whole in- 
terval of the cutoff parameter G means an essentially 
non-Fermi-liquid (NFL) picture. A similar mechanism 
of NFL behavior in magnetic Kondo lattices was pro- 
posed in Refs0,0 where a soft boson was obtained with 
the characteristic energy tU(C) = \G\. Note also that a 
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FIG. 1: Scaling trajectories -Ef{^) for V{0) = 0.1, G(0) = 
0.1, N — 2 in the incoherent case as compared to the 
curve \C\/D — exp(— 5) (dashed line). The parameter 
values (for the curves from below to above) are Ef{0) — 
-0.08,-0.1,-0.14. 




FIG. 2: Dependences V{^) for the same parameter values as 
in Fig.l (for the curves from below to above) 



regime with the rate of the order parameter fluctuations 
^/T(I>{T) oc T is typical near a quantum phase transi- 
tion |2l| . Our situation is reminiscent of this regime in 
the sense that the characteristic valence fluctuation fre- 
quency ma.x{pV^{T), \Ef(T)\) is proportional to T (after 
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FIG. 3; Scaling trajectories -£/(0 for V{0) = 0.05,0(0) = 
0.05, iV = 6,Ef{0) = -0.09,-0.11,-0.15 in the coherent 
case. 
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FIG. 4: Dependences V{£,) for the same parameter values as 
in Fig. 3 (for the curves from below to above) 



a natural replacement |C| — > T). 

Physically, the regime where a typical energy scale 
is just the temperature means a classical (Maxwell- 
Boltzmann) electron gas (interacting with local mo- 
menta): heat capacity is approximately constant etc. Of 
course, in a pure form this behavior is never observed 
since, according to our results, it takes place only in a re- 
stricted temperature interval. However, we have a strong 
deviation from a simple scaling picture where we just 
enter strong-coupling regime, a characteristic "Kondo" 
temperature being T-independent. In particular, the 
Wilson number is nearly constant, but differs consider- 
ably from that in the singlet state (Ig/ — — oo). 

A similar "marginal" region in the dependence Ef(C) 
occurs not only for the Kondo, but also for IV state near 
the critical line (in such a situation, the dependence w{C) 
has a shallow minimum, and V{C) a sharp maximum, see 
Figs. 1-2). 

To estimate the Kondo temperature we can use Hal- 
dane's arguments for the Anderson model with N = 2 
|l3l | . The generalization to arbitrary TV can be performed 
as (cf. Ref. 13) 
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This expression is formally based on perturbation theory 
(two-loop scaling). However, Haldane noted that replac- 
ing in this formula both D and Ef by the characteristic 
energy scale T* = — C, which is determined from the 
equation C = Ef{C), yields the correct estimation for 
the Kondo temperature. As demonstrated above, in the 
presence of exciton effects (Falicov-Kimball interaction 
G), a situation is possible where this equation holds ap- 
proximately in a whole energy interval. The energy scale 
T* where the "marginal" regime starts is considerably 
changed by the exciton effects. 

In the coherent case the last term in Eq. Q results in a 
smearing of the singularity, especially for small N. How- 
ever, with increasing N the dependence Ef{C) (Fig. 3) 
becomes qualitatively similar to that in the incoherent 
case. On the other hand, the dependence V{C) (Fig. 4) 
is essentially modified even for — > oo. At C = Ef{C) 
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we have 

1 = -{N-l)pV^{C)/C-n[G{Q)/V{0)]dV{C)/dC (14) 

Solution to this Riccati equation is obtained in terms 
of the imaginary-argument Bessel and Macdonald func- 
tions Ip{x) and Kp{\{N - l)pC\'^''^) with x = 2\{N ~ 
l)pV'^{0)C\^/^/[nG{0)], p = 0, 1. For large N, the quan- 
tity dV{C)/dC and consequently V^{^)/Ef{^) turn out 
to be practically ^-linear near the maximum of w{C). 
Thus we have a "classical" electron liquid with singular 
interactions which have logarithmic energy dependences. 
With further decreasing |C| we have from Eq.© 

dV{C)/dC oc \C\-^/^ = \Ef{C)\-^/^ (15) 
so that the correction to V'^{C)/Ef{C) is proportional 

to |C|-l/2. 

The d — f Coulomb interaction can strongly renormal- 
ize the hybridization, which leads to the increase of the 
characteristic energy scale. In the incoherent regime it is 
a width of the resonance, in the coherent one the width 
of (indirect) gap or pseudogap The renormalization 
of the fluctuation rate pV^{C) can be very strong (about 
order of magnitude for realistic parameters). The corre- 
sponding temperature dependences can be found in both 
the regimes by the RG approach with the replacement 
\C\^T . 

In standard treatments of HF systems, one picks 
usually anomalous magnetic contributions to thermo- 
dynamic properties and compares them with exact re- 
sults in the one-impurity Kondo problem. In partic- 
ular, the dependences of the crystal-field level width 



from both ineleastic neutron scattering and nuclear mag- 
netic resonance (NMR) have the form r(T) oc T^/^. In 
the Kondo resonance model such a behavior takes place 
above the Kondo temperature Tk, but in some cases 
(CeBe, CePdsBo.e, YbBcia) the dependence T^/^ ^^y-^g 
place at very low temperatures of a few K 22]. Further, 
the characteristic energy scale from linear heat capacity 
term 7 is of order of tens of Kelvins, whereas the tem- 
perature where 7 starts to deviate from constant is just a 
few of Kelvins. Thus there exists no unique energy scale. 
We have demonstrated that, indeed, the infrared behav- 
ior can be essentially different from the simple Anderson 
model owing to spin dynamics (see RefsiTi il^) and charge 
fluctuations (exciton effects) considered in the present 
work. 



Of course, the estimations performed are qualitative 
since they are based on the continuation of perturbative 
Gellmann-Low scaling function to the strong coupling re- 
gion. At the same time, the statement that exciton effects 
cannot be described by universal temperature indepen- 
dent Tk seems to be reliable itself. Recently, direct ways 
of observing the Kondo resonance (STM) were proposed 
[23l |. As we know, they have been not yet applied to 
the IV systems. It would be interesting to compare the 
results for Tk of these new experimental methods with 
those from investigating thermodynamic properties. 
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